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It is pointed out recentlyEI that the v = 1/m quantum Hall states in bilayer systems behave like 
easy plane quantum ferromagnets. We study the magnetotransport of these systems using their 
"ferromagnetic" properties and a novel spin-charge relation of their excitations. The general trans- 
port is a combination of the ususal Hall transport and a time dependent transport with quantized 
time average. The latter is due to a phase slippage process in spacetime and is characterized by 
two topological constants. 
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Recent technological advances in producing high mo- 
bility bilayers have provided opportunities for studying 
quantum Hall (QH) systems with internal degrees of free- 
dom ||. In this case, the internal degrees of freedom is 
a pseudo-spin 1/2 labelling the upper and lower layer. 
There are experimental evidence g indicating that un- 
der appropriate conditions, QH states (at v = 1/2 and 
v = 1) can be stabilized by correlations between layers 
and have little to do with tunneling between them. It is 
conceivable that similar correlated states will be found 
at lower filling factors. Bilayer QH states with v = 1/m, 
(m odd) are believed to be the so-called (mmm) state. It 
is realized recently that they are very unusual QH fluids. 
Wen and Zee Q showed that these states can produce 
Josephson like tunneling between the layers. Very re- 
cently, a number of authors have pointed out that 
these states behave like easy plane quantum ferromag- 
nets and have studied possible phase transitions in these 
systems. 

The purpose of this paper is to study the magneto- 
transport of QH ferromagnets (QHFs). The difficulty of 
the study arises from the fact that the usual procedure of 
extracting conductivity tensor (say, using flux insertion 
arguments or Chern-Simon type of theories) breaks down 
when applied to the (mmm) state Here, we take a 
new approach which amounts to studying the hydrody- 
namics of QHFs. Our approach recovers the known re- 
sults for ordinary QH fluids, and reveals a new transport 
mode in QHF. A general QHF transport is a combination 
of usual QH trasnport and a genuinely time dependent 
transport with quantized time average. The latter corre- 
sponds to a phase slippage process which turns out to be 
a continuous nucleation of coreless vortices in spacetime. 
The conducitivity tensor for the time average currents is 
characterized by two topological constants specifying the 
nucleation process. It is interesting to note that there 
is a close analogy between QHF and superfluid 3 He-A. 
They have similar order parameters, topological excita- 
tions, and vortex nucleation processes. 

To determine the magnetotransport of a bilayer sys- 
tem, we need to find the current densities J M for given 



electric fields E M in each layer, where [i =f and J, denote 
the upper and lower layer respectively. A general con- 
figuration {E^} is a sum of identical and opposite fields 
(E T ,E i ) = E+(l,l) + E_(l,-l);E ± = ±(E T ±E X ). The 
cases (E+ / 0, E = 0) and (E + = 0, E_ ^ 0) will be 
referred to as the "in-phase" and "out-of-phase" mode. 
Their corresponding current responses will be denoted by 
a subscript "+" and "— " respectively. In linear response, 
it is sufficient to consider these modes separately. The 
"in-phase" mode is the usual QH transport, with current 
response 

(J,|Jj) + = ^BxE + , (J T -J x ) + = (1) 

where B is the direction of the magnetic field, taken to 
be normal to the layers in our discussions. The "out-of- 
phase" mode is unique to bilayer systems. For QHFs, we 
shall show that d.c. fields E_ will generate time depen- 
dent currents. The sum (Jj + Jj.)_ oscillates with period 
T = 2T a , where T a is the Josephson period eAV/h, and 
AV is the voltage that generates E_ . The factor of 2 is a 
reflection of the double- valueness of the pseudo-spin 1/2. 
The time average current response (denoted by a "bar" ) 
can be determined from general argument and are given 
by 
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(J T + J;)_ = nAV—B x E_, (J T - J x )_ = B x 

(2) 

where ua and jib are integers. By studying a number of 
examples consistent with the driving force E_, we find 
UA = 0, riB = 1. The sum of eq.(|l]) and (||) then gives 

These examples also reveal an interesting oscillatory be- 
havior of the currents. The time average of (Jj + JjJ- 
over the half period T Q is ±i/(e 2 //i)B x E_ , with alternat- 
ing sign every half period. The time average of (Jj — Jj.)- 
is (v/2)(e 2 /h)B x E for all half periods. 
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To derive these results, we recall that the Hamiltonian 
of the bilayer system in the weak tunneling limit is @ 



H 



2m J 



I A) ip — Ami(r) 



+ 2 / [p( r )p( r X r - r') + m 3 (r)m 3 (r>i(r - r')] (4) 

where r = (x,y), VVj(r) creates an electron at the p layer 
at r, rhi = hTrip + cr^ , (i — 1, 2, 3), p = Trr(j + 4> and both 
v, Vi > 0. The pseudo-spin axes Xi,X2,X3 are chosen so 
that the symmetric and antisymmetric state (with energy 
difference A) are represented by Q) and (-t). They have 
no relations with the real space directions x, y, z. 

We shall consider the (mmm) state in the correlation 
limit A << e 2 /e £ ~ e 2 /e Q d, where d is the spacing be- 
tween layers, and e Q is the dielectric constant. We be- 
gin by expressing their wavefunctions (in the symmetric 
gauge) in spinor form 

JV N 

2—1 i>j 

where z = x + iy and 

c = o = ici (::^ 2 ) e - x/2 (6) 

is a spinor along £ = X3Cosf? + sin#(xiCOS</> + X2sin0). 
Since v\ > 0, it is energetically more favorable for I to lie 
in the X1-X2 plane. A will line up £ along xi. However, 
in the correlation limit, it remains small compared to 
the correction energy. The symmetric and antisymmet- 
ric states can therefore be regarded as essentially degen- 
erate a view that we shall adopt in our subsequent 

discussions for simplicity j8j. 

If uniform potentials V\ and Vj are applied to in the 
upper and lower layer, u and v will evolve as e leV r*/ ri and 



a ieVit/h 



(e > 0), meaning 



(7) 



Eq. (0) implies means that £ precesses about X3 with an- 
gular frequency eV/h. This immediately implies an "un- 
stable" situation if V is different at two points far apart, 
(i.e. E_ ^ 0). The I vectors at these points will assume 
their ground state configurations, (lying in X1-X2 plane) 
but rotate about x 3 at different rates. As a result, more 
and more winding in the spin texture is produced which 
can only be relieved by phase slippage processes. As we 
shall see, during phase slippage, I will be pulled away 
from the X1-X2 plane every now and then in some region 
in space. While this costs nucleation energy because of 
the last term in eq.(f|), it can not be stopped because its 
energy cost will evenetually be surpassed by the gradient 
energy generated by E_. 



Nonuniform textures implies spatially varying spinors. 
However, to stay in the lowest Landau, Q must be analytic 
in z, (or a function of d z ). For simplicity, we shall from 
now on focus on the "quasihole" case where ((z) is an 
analytic function of z. In general, £ can be expanded 
in powers of z. The simplest example of a nonuniform 
analytic spinor is £(z) = a + j3z, where £ reduces to 
spinors a and (3 as z approaches and 00. It also reduces 
to the usual quasihole excitations if a oc f3. For arbitrary 
C, eq.(||) implies that 



2m(r) = p(t)£(t) 



(8) 



This is a statement of "full" ferromagnetism and is im- 
portant for our subsequent discussions. 

The relations between density and spin fluctation can 
be obtained from the "plasma analog" B, where one 
interprets 1 |*([r, /j])\ 2 as the partition function of 
a classical 2D plasma with charge m in a background 
potential V b (z) = -\z\ 2 /2 + ln|C| 2 . This implies that 
the electron density is p(r) = = — (47tto) _1 V 2 VJ,. = 
(27rm) _1 — (47rm)~ 1 V 2 ln|('| 2 . Because of the analyticity 
of C, the density change can be expressed in terms of a 
deformation field @ u = -<+VC/|C| 2 + c - c - 



P(r) = 



1 



1 



27TTO Airm 



-z • V x u 



(9) 



If follows from eq.(|) that Q fjj 

z • V x u = -(z x V) • Vx + ie Q/ 3 7 ^ Q (z x V)lp ■ V£ 7 . 

(10) 

The first term vanishes unless x is a vortex, in which case 
it is a (5-function Q . 

The total electric current J = Jj + J| can be obtained 
by differentiating the number density eq.(|^) and using 
the continuity equation (— e)p = — V • J, 



1 e 2 e 

J a = — — (z x V) a U - 

m h vktyi 



£-{zxv) a exd t i (11) 



where we have made use eq.(Q). The first term is the 
usual QH transport, eq.(|l|). (Note that E + = — Vf7 and 
that eq.(||) is derived for B = — z). The second term is 
due to the motion of the spin texture, which is generated 
by E_. To simplify discussions and formulas, we shall 
from now on focus on the out-of-phase mode, i.e. VJ7 = 
0, W 7^ 0. We shall then drop the first terms in eq.(|Tl|) 
and eq.®. 

Consider now the time average of the current in y over 
a period T, (7f + I^) y = T^ 1 L dt J J y dx. The integral 
is evaluated for an arbitrarily y coordinate, with upper 
and lower limit x a and x b denoting the edges of the sam- 
ple. From eq.(ll) we have 



(M+h)y 



1 



47T77J T 



dt I dxl- dj x dj 



(12) 
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The quantity in the bracket in eq.([12|) (denoted as A) is 
the area on an unit sphere 62 covered by the "spactime 
rectangle" T with area T x L, L = (xb — x a ). A non- 
vanishing A means nucleation of vorticity in spacetime 
[fit(Vxu)^0]. 




alone cannot determine I<\ y and I 
|), and (|ic|), we have 



Iv 



2m(r) 



1 



-I 
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d x £ x d v £ 



2<9 t m = 



* 9 t £ 

znm 4nm 



— V Q f (z x V)J x d t 



From 



(13) 



(14) 



Eq.(J3) says that spin fluctuations will generate a 
pseudo-spin current 



Kq, = (z x V) i x c\i 
47rm 



(15) 



which is related to the charge currents as (— e)x3 • K Q = 
J] a — J [a- The time average of current difference in y is 



-e 1 

47T771 T 



dt 



dx X3 • d x £ x dti 



(16) 



The quantity in the bracket (denoted as B) is the area 
in 52 covered by T projected onto the X1-X2 plane. 

To evaluate A and B, we consider a rectangular sam- 
ple (see fig.l) where the potentials (V\, VjJ at x a and Xb 
are (0,0) and (V/2, — V/2) respectively. This is the "un- 
stable" situation we discussed before, where I precesses 
about X3 with angular frequency eV/Ti at Xb and remains 
stationary at x a p4[ . In the steady state, the textural 
motion must be periodic as E_ forces a constant rotation 
on ( at Xb- Note that £ at Xb returns to itself after time 
T = h/eV, whereas the period of £ (and hence the en- 
tire texture £(x,y,t)) is T = 2T D . The factor of 2 is due 
to the double valueness of £. An example of the steady 
state evolution of £ is shown in the spacetime plot (x-t) 
in figure 1. The evolution is a repetition of the pattern 
inside the rectangle T x L, (abb" a"). Since the bound- 
ary of the rectangle maps onto the equator of £2 twice, 
we then have A = —AirriA; B = —2irriB, where tia and 
ub are integers depending on how S2 is covered. Eq.jl^) 

n A (l/m)(e 2 /h)(V/2), 



and ( |16| ) then gives (Ij + Ii) y 

(J T - Jjj = n B (2m)- 1 (e 2 //i)(y/2), which is eq.(|) since 
E_ = -|V(Vt - Vi) = -^Vy and B = -z. 

To determine the topological integers n^, we con- 
struct explicit examples for the out-of-phase current re- 
sponse. Since we are dealing with rectangular geometries, 
it is more convenient to use the Landau gauge. The QHF 
now reads [051 



ncK)ii(^-^) r 



where w = e az , a is the ratio of the magnetic length to 
the sample length in y, which can be taken as 1 with- 
out loss of generality. To match the Josephson relations 
eq.(0) at the boundaries x a and Xb, (now taken to be —00 
and +00), we choose a spinor of the form 



CM = A 



1 



JeVt/2h 
-ieVt/2h 



(18) 



where A is a constant setting the scale of textural destor- 
tion in real space |]l6| . Its textural evolution in the space- 
time plane (x-t) (for a given y) is shown in fig.l. £ be- 
comes uniform every half period T a = h/eV. The topo- 
logical areas A and B for each half peirod can be calcu- 
lated from eq.(18|) and are independent of A. A is found 
to be =f27t, with alternating signs every half period. B is 
found to be — it for all half periods. This means over the 
full period T = 2T a , A = and B = —2ir, corresponding 
to ua — 0,n.B = 1. We have then derived eq.(||) from 
this specific example. I have also examined a number 
of other spinors satisfying the same Josephson boundary 
condition at infinities. Their time average responses are 
identical to that of eq.(|l8|). I therefore think that these 
values of ua and fig have general validity, despite the 
lack of a general proof which is certainly desirable. 
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Figure Caption 

Figure 1 : The I vectors at r is represented by a line 
emerging from r. An solid (empty) square is placed at r 
to indicate that I has a positive (negative) 23 component. 
The boundary of the rectangle T — (abb' a') maps onto 
the equator T on the unit sphere S2. The lines cc', dd', 
ee', c'c", d'd", and e'e" map onto loops C, D, E, C, D, E. 
This covering gives A — 2tt, B = —tt for the lower rectan- 
gle abb' a' and A = —2ir, B = —tt for the upper rectangle 
a'b'b"a". If the loops C, D and E are traversed in direc- 
tions opposite to those depicted, then A = 2ir,B = n for 
a'b'b" a" . The textural pattern is calculated from eq.jT^) 
with T = l,y = 2, A = 8. 
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